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Cryptochromes are flavoproteins that drive diverse developmental light-responses in plants and participate 
in the circadian clock in animals. Plant cryptochromes have found application as photoswitches in 
optogenetics. We have studied effects of pH and ATP on the functionally relevant photoreduction of the 
oxidized FAD cofactor to the semi-reduced FADH radical in isolated Arabidopsis cryptochrome 1 by 
transient absorption spectroscopy on nanosecond to millisecond timescales. In the absence of ATP, the yield 
of light-induced radicals strongly decreased with increasing pH from 6.5 to 8.5. With ATP present, these 
yields were significantly higher and virtually pH -independent up to pH 9. Analysis of our data in light of the 
crystallographic structure suggests that ATP-binding shifts the pKa of aspartic acid D396, the putative 
proton donor to FAD*", from —7.4 to >9, and favours a reaction pathway yielding long-lived aspartate 
D396". Its negative charge could trigger conformational changes necessary for signal transduction. 

Cryptochromes (CRYs) exhibit diverse biological functions, yet share an essential, non-covalently bound 
flavin adenine dinucleotide (FAD) cofactor^'^. In mammalian CRYs, where the FAD-binding pocket can 
be open and dynamic^'^, replacement of FAD by a carbazole derivative affected the length of the circadian 
period^'^, suggesting therapeutic potential. The active redox states of FAD in CRYs are still a matter of debate and 
likely differ among organisms and/or for the given protein function^"^. 

In Arabidopsis thaliana cryptochrome 1 (AtCRYl), the semi-reduced (semiquinone) neutral FAD radical 
(FADH ) appears to be the signalling state^. This signalling state is formed by photoreduction and subsequent 
protonation of the resting state, the fully- oxidized FAD (FADox), which strongly absorbs blue light (£445 = 
11 300 M"^cm"^)^°. By analogy to the structurally related DNA repair enzyme photolyase^\ in which fully 
reduced FAD is replenished by so-called photoactivation, the reducing electron in AtCRYl is transferred to 
the excited FAD (*FADox) through a cascade of tryptophan residues'^''' (W400, W377 and W324, see Fig. 1; we 
denote them by TrpiH, Trp2H and Trp3H, respectively, when the protonation state is relevant) and possibly also a 
tyrosine (TyrOH)^^ residue. The quantum yield of FADH" Trp' radical pair formation in isolated AtCRYl was 
estimated to be as low as —2%^^. The facile reduction of Trp3' and TyrO' by extrinsic reductants stabilizes FADH' 
by preventing recombination of the radical pairs. The mechanism by which the FADox FADH' transition 
might be sensed by signalling partners remains to be established. 

Both plant and animal CRYs have been reported to bind ATP and exhibit autokinase activity in vitro^"^'^^, 
though the implications of such enzymatic activity have yet to be elucidated in vivo. Furthermore, a prolonged 
lifetime of FADH', the putative active redox form of FAD, and a higher resistance to denaturation and proteolysis 
have been reported for plant cryptochromes in the presence of ATP^^'^^. Isolated AtCRYl binds ATP with a 
stoichiometry of 0.9-1.0 ATP/protein^^'^° and a dissociation constant on the order of 10~^ ]y[i4,i6,2o -pj^^ photo- 
lyase homology region (PHR) of AtCRYl, lacking the poorly structured C-terminal region (CCT), was 
CO -crystallized with a non-hydrolysable ATP analogue (adenylyl-imidodiphosphate, AMP-PNP). The crystal 
structure^° revealed that the ATP-binding site is located near the FAD cofactor. CRYs show a strikingly similar 
structural fold to photolyases^\ consisting of N-terminal a/p and C-terminal a-barrel domains. FAD is buried in 
the a-barrel domain. Both CRYs and photolyases have a pocket leading to the FAD cofactor; ATP binds in the 
pocket of plant CRYs, whereas photolyases recognize pyrimidine-dimer lesions and restore normal bases at this 
site. When bound to AfCRYl, the adenine moiety of ATP is tightly sandwiched between leucine L296 and tyrosine 
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Figure 1 | Tertiary structure of AtCRYl. (a) Overall structure of the photolyase homology region in Arabidopsis thaliana cryptochrome 1 with a bound 
non-hydrolyzable ATP analogue, adenylyl-imidodiphosphate (AMP-PNP), adapted from the crystal structure by Brautigam et al?^. The N-terminal a/p 
domain exhibits a Rossman fold; the FAD cofactor (yellow) and AMP-PNP (pink) are bound to the a-helical domain. The carboxy-terminal CCT domain 
comprising another —220 amino acids in full-length AfCRYl is absent in the crystallized protein^°. (b) Detailed view of the surroundings of the FAD 
cofactor. The flavin-reducing tryptophan triad W400 (TrpiH), W377 (Trp2H) and W324 (TrpsH)^^ is shown, together with the putative proton donor 
aspartic acid D396, and key residues interacting with the adenine moiety of AMP-PNP: leucine L296 and tyrosine Y402 (note that the al6-helix 
interconnects Y402 with D396 and W400). The H-bond between the aspartic acid D396 and the backbone carbonyl oxygen of methionine M381 
interconnecting the helices al5 and al6 and the ultrafast proton transfer pathway^^ between W400 and N5 of FAD through D396 are indicated by 
magenta dashed lines. The electron transfer (ET) steps are shown by red arrows. 



Y402, which are highly conserved among plant CRYs. ATP phos- 
phates are anchored by positively- charged residues on the protein 
surface. 

Here we demonstrate by transient absorption spectroscopy that, in 
addition to the aforementioned effects, ATP binding to AtCRYl 
leads also to an important increase in the yield of light-induced 
FAD'"-TrpH"^ and FADH -Trp' radical pairs, and, accordingly, 
also in the yield of the metastable isolated FADH' radical. In the 
absence of ATP, we found that the efficiency of CRY photoreduction 
decreases substantially with increasing basicity between pH 6.5 and 
8.5. ATP binding remediates this decrease, presumably by shifting 
the pH- sensitivity region towards the basic (to pH > 9.0). 

Detailed analysis of the transient absorption kinetics observed in 
our study shows that both pH and ATP binding affect the protona- 
tion of FAD ", formed by photo-induced electron transfer (ET) to 
*FADox from W400. Our data suggest that the observed pH-effect 
reflects the protonation state of D396(H), which is positioned oppos- 
ite the N5 of the isoalloxazine ring of the FAD, and that ATP binding 
up-shifts the piCa of this aspartic acid. D396(H) is widely accepted to 
be the proton donor in the formation of FADH' from FAD " in 
AtCRYl however, there is currently a controversy on whether 
FAD " protonation occurs in an ultrafast manner (<1 ns)^^ or on a 
microsecond scale^^ in plant crypt ochromes. Our results indicate that 
both pathways exist and are affected by pH and by ATP. 

Results 

Transient absorption studies. Addressing the question of whether 
or not ATP directly enhances the photoreduction of FADqx in CRYs 
(in addition to the reported stabilization of the FADH" form^^'^^), we 
have monitored the yields and kinetics of species involved in FADox 
photoreduction by transient absorption spectroscopy in the presence 
and absence of ATP. 

Based on the absorption spectra of the expected transient species 
(Fig. 2), we have chosen three basic spectral regions to monitor the 
photoreactions subsequent to FADox excitation: 400-410 nm. 



—450 nm and —560 nm. At —450 nm, FADqx absorbs much stron- 
ger than both FAD " and FADH" and the absorption of other 
expected species is negligible. Hence the major processes observed 
at this wavelength are FAD" "/FADH" formation (initial bleaching) 
and their recombination with tryptophan radicals back to FADox 
(recovery of bleaching). The region around 560 nm allows monitor- 
ing of TrpH "^ deprotonation to Trp" (decrease of the initially formed 
absorption) and the protonation of FAD " to FADH" (absorption 
increase). Finally, at 400-410 nm, formation of FAD" is accompan- 
ied by an absorption increase, and the subsequent protonation of 
FAD " to FADH" is reflected by decay down to a slightly negative 
absorption difference (AA) signal (FADH" absorbs less than FADox 
in this region). 

In order to obtain the best possible signal-to-noise ratios of some- 
times weak transient absorption signals, we have used truncated 
AtCRYl that lacks the flexible C-terminus and corresponds to the 
N-terminal PHR. This protein proved to be more stable and samples 
of higher concentration, higher purity and better batch-to-batch 
consistency could be produced in sufficient amounts. We have con- 
ducted several control experiments with full-length (wild-type) 
AtCRYl to verify that the photochemistry of the two proteins is 
not substantially different (Supplementary Results, Supplementary 
Fig. 2). Indeed, we have not found any significant differences in their 
photoreactions and thus we conclude that the C-terminus has vir- 
tually no effect on the primary photophysical and photochemical 
events in AtCRYl. 

The time resolution of the transient absorption set-up was limited 
to —10 ns in order to obtain a good signal-to-noise ratio. Quantum 
yields of transient radical pairs were determined using the 
[Ru(bpy)3]^^ actinometer^^. 

All kinetic traces in Fig. 2 clearly show that ATP substantially 
increases the yield of the photogenerated transient species. In addi- 
tion, ATP changes the evolution of their relative concentrations in 
time (best seen at 562 nm, but also at 457 nm), which is a clear 
evidence of a mechanistic change in the course of the photoreaction. 
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Figure 2 | Nanosecond-microsecond kinetic traces for AfCRYl PHR. Transient absorption changes were recorded at 457, 562 and 403 nm upon 
excitation of FAD ox by a picosecond laser pulse (3.0 mj at 355 nm, at time zero) in the absence (black traces) and in the presence of 1 mM ATP (red 
traces). The poorer signal to noise ratio at 403 nm is due to an inherently worse stability of the monitoring light source at this wavelength. The samples 
containing 47 fiM protein, 0.05 M phosphate buffer of pH 7.4 and 0.5 M NaCl were kept at 1-2 °C. Top-right segment: absorption spectra of FAD in three 
of its redox forms (FADox> FAD" and FADH) and other species potentially contributing to the absorption changes. The FAD ox spectrum is that of our 
AfCRYl PHR, normalized to £445 = 11 300 M"^cm"^ (from Ref. 10). The FADH' spectrum has been extracted from a spectrum of a mixture of FAD ox and 
FADH' of the same protein (see Supplementary Fig. 7 for more details). The FAD'~ spectrum is that of FAD'" in a mosquito cryptochrome (A^RYl) 
adopted from Liu etaV^. TrpH'^ and Trp' spectra are adopted from Solar etaV^, and TyrO' from Giese et aV^, normalized to e {Irnax) — 2 600 M~^cm~^ 
(from Dudley Bryant et aV^). 



Two major kinetic phases are apparent in all signals (see fit in 
Supplementary Fig. 8a): The faster phase has a time constant t ~ 
200 ns and is assigned to deprotonation of TrpH"^ (contributing 
most significantly at 562 nm) parallel to recombination of FAD " 
and TrpH"^ (contributing most significantly at 457 nm). The slower 
phase (t ~ 1.5 |is) is assigned to the protonation of FAD " to FADH" 
(absorption increase at 562 nm and decrease at 403 nm). Both 
phases are relatively more pronounced in the presence of ATP than 
in its absence. The absorption changes at the end of the time scale of 
Fig. 2 (10 |is) are assigned to mostly FADH' - Trp" radical pairs 
(—80%) with a small admixture of FAD " - Trp' pairs (—20%; see 
Supplementary Information, Section 8). 

Signals from ATP-free samples in Fig. 2 showing only a little 
evolution in the given time window seem to be consistent with the 
recently theoretically predicted^^ ultrafast protonation of FAD " by 
D396(H) and a subsequent protonation of D396" by TrpjH"^ {i.e., 
formation of an FADH" - Trp" pair already on the picosecond time- 
scale). On the other hand, in samples containing ATP (and to a 
smaller but non-negligible extent also in samples lacking it) we 
clearly observed a much slower FAD " protonation kinetics (t — 
1.5 jis; close to the 1.7 |is value observed in experiments on algal 



cryptochrome^^) preceded by TrpH"^ deprotonation with an appar- 
ent T — 200 ns (comparable to deprotonation of the terminal TrpH" ^ 
in DNA photolyase from E. colf^ indicating a reaction pathway 
different from the one predicted by the theoretical study^^). 

Striving to better resolve and understand these protonation/ 
deprotonation processes, we have recorded kinetic traces at 
562 nm at different pH values (Fig. 3). The signals revealed another 
remarkable phenomenon: the yield of photogenerated species is very 
sensitive to pH change in the absence of ATP, diminishing dramat- 
ically with increasing basicity from pH 7.0 to 8.0. The signals from 
samples containing ATP, however, remain virtually unchanged 
within this pH range and decrease only at the highest pH tested 
(9.6). The pH-dependence of signal amplitudes from ATP-free sam- 
ples after 10 |is (with major contribution to AA from the FADH" - 
Trp" radical pair) can be fitted by a titration curve for a monoprotic 
acid with piCa — 7.4 (Supplementary Fig. 3a). 

The evolution of the transient species on time scales longer than 
10 |is was monitored on a set-up with lower time resolution but 
higher amplitude sensitivity than the set-up used for Figs. 2 and 3. 
At all wavelengths studied (410, 450 and 560 nm), most of the 
detected absorption change decayed in a few milhseconds (Fig. 4). 
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Figure 3 | Nanosecond-microsecond kinetic traces for AfCRYl PHR at different pH values. Transient absorption changes were recorded at 562 nm 
upon excitation of FAD by a picosecond laser pulse (3.5 mj at 355 nm, at time zero) in the absence (black traces) and in the presence of ATP (red traces). 
While 1 mM ATP was largely sufficient to saturate the ATP effect on the signal amplitude at pH < 8, more ATP had to be added to the sample at pH 9.6 in 
order to reach full saturation (possibly due to the fact that ATP hydrolyses rapidly at high pH). The samples contained 45 jjM protein, 0.05 M phosphate 
buffer of pH 7.0 to 8.0 or 0.05 M Tris buffer of pH 9.6 and 0.5 M NaCl and were kept at 1-2°C. 



We attribute this decay essentially to recombination of the predom- 
inating pair of neutral radicals FADH" - Trp". This recombination 
seems to be accelerated by ATP (from t = 8.0 ms to 2.2 ms; see 
Supplementary Information, Section 4 for possible origins of this 
effect). The presence of even longer-lived absorption changes with 
spectral features reminiscent of FADH" indicates that part of FADH' 
was stabilized by reduction of Trp" by an intrinsic reductant different 
from FADH", most likely a neighbouring tyrosine residue (see Fig. 4). 
FADH" - Trp" recombination can be virtually completely suppressed 
by addition of an extrinsic reductant (Supplementary Fig. 5). As a 
result of efficient reduction of Trp", the overall quantum yield of the 
isolated metastable flavin radicals (FADH"/FAD"") can attain —14% 
in the presence of ATP (Supplementary Information, Sections 4 
and 8). 

Change in the UV/Vis spectrum associated with ATP binding and 

pH. While checking the quality and the state of our samples by UV- 
Vis spectroscopy, we have noticed that binding of ATP by A^CRYl is 
accompanied by a reproducible change in the UV-Vis spectrum of 
FADox (Fig. 5a and Supplementary Fig. 2b). This effect is similar to 
the reported spectral change accompanying the substrate binding by 
a DNA photolyase containing an oxidized FAD cofactor^^'^°. While 
the molar absorption coefficients (e) around the maximum in the 
blue spectral region (445 nm) slightly increase in the presence of 
ATP (by -400 M"^cm"^ at pH 7.4) and the maxima exhibit a 



slight red shift (—1-2 nm), the two maxima in the near-UV (361 
and 371 nm) shift towards shorter wavelengths (by —4 and —1 nm, 
respectively), and their amplitude ratio reverses. The ATP 
dependence of the spectral change can be described by a 
dissociation constant iC^ = 1.4 ± 0.6 jiM (see Supplementary Fig. 4). 

Interestingly, we have observed a similar spectral change when 
lowering the pH in the absence of ATP (Fig. 5b). Samples at lower 
pH showed blue shift of the UV bands with respect to samples at 
higher pH, but the red shift of the blue band was not so distinct. 
Again, like in the case of kinetic traces, also this pH-effect faded away 
upon ATP addition (in the studied pH range, i.e., between pH 6.5 and 
9.6) and there was no significant difference in the spectra of samples 
saturated with ATP at different pH values. The spectral shift with pH 
observed for ATP-free samples is a strong indication of a change in 
the protonation state of an amino acid in the vicinity of the FADqx 
cofactor. Plotting the absorption difference at 386 nm (minimum of 
the negative peak in the difference spectrum, inset Fig. 5b) against pH 
(Supplementary Fig. 3b) yields a titration curve of a monoprotic acid 
with the same piC^ as the one obtained from the 562 nm- signal 
amplitudes after 10 |is, i.e., 7.4 (Supplementary Fig. 3a). 

Discussion 

Three major previously unreported phenomena were discovered in 
the present study of flavin photoreduction in isolated plant CRY on 
nanosecond to millisecond timescales: 
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Figure 4 | Millisecond kinetic traces recorded for AfCRYl PHR. 

Transient absorption changes upon excitation of FAD ox by a picosecond 
laser pulse (2.2 mj at 355 nm, at time zero) were recorded at 410, 450 and 
560 nm in the absence/presence of 1 mM ATP. The samples containing 
35 fiM protein, 0.05 M phosphate buffer of pH 7.4 and 0.5 M NaCl were 
kept at 1-2 °C. Note that despite ATP -accelerated recombination, the 
signal amplitudes remaining after 100 ms show still —70% higher amounts 
of FADH' in samples containing ATP. At 410 nm, the signal bleaching due 
to the difference in e of FAD ox and FADH' (see spectra in Fig. 2) is largely 
compensated by the absorption of the newly formed TyrO' radical. 

(i) In the absence of ATP, the yield of transient species formed 
upon photoexcitation of FADox dramatically decreases with 
increasing pH in the accessible pH range (6.5-9.6). The yield 
evolution can be fitted with a titration curve for a monoprotic 
acid of piCa ~ ^-4. 

(ii) Addition of ATP strongly increases the yield of transient spe- 
cies and makes it pH -independent up to pH ~ 9. Moreover, 
ATP induces substantial qualitative changes in the species kin- 
etics. 

(iii) In the absence of ATP, the near-UV absorption spectrum of 
FADox in AtCRYl changes with pH and this evolution also fits 
a pKa ~ 7.4. Addition of ATP and lowering the pH impacts 
the spectra in a similar way. When saturated with ATP, the 



spectrum of FADqx in AtCRYl no longer changes with pH (in 
the studied pH range, which was limited by fast aggregation of 
the protein at pH < 7 and by ATP hydrolysis at pH > 9). 

In the following, we put our results into context with the current 
literature and propose a mechanistic model for FAD photoreduction 
combining photolyase-like ET through the tryptophan cascade^ ^ 
with the theoretically predicted^^ D396(H) -assisted ultrafast PTs 
from TrpiH + to FAD " (Fig. 6). 

Upon excitation of FADox, ultrafast (0.4 ps^^) ET from the nearby 
TrpiH forms the first charge separated radical pair FAD " - 
TrpiH"^. Subsequent ET through the tryptophan cascade (presum- 
ably in <1 ns) separates the charges farther. The resulting pair 
FAD " - Trp3H"^ is stabilized by tryptophan deprotonation to the 
solvent within —200 ns (observed in the present study), disfavouring 
further unproductive back ET in the FAD " - TrpH"^ pairs. So far, 
the reaction scheme is analogous to photoreduction of FADH" to 
FADH" in DNA photolyase^\ However, an additional PT to 
FAD " occurs in plant CRY to yield the putative signalling state 
FADH'. Based on structural considerations, FTIR data and transient 
absorption data^^"^^ it appears estabhshed that FAD " is protonated 
by a nearby aspartic acid (D396 in AtCRYl) in ~ 1.5 |is (observed in 
our study; 1.7 |is was reported for CRY from Chlamydomonas^'^). 

In CRY/PL proteins that contain asparagine in the position cor- 
responding to D396 in AtCRYl (or even in the D396N mutant of 
AtCRYP^), FAD " can get eventually protonated via alternative 
pathways^^. Nevertheless, the only time-resolved study of FAD " 
protonation in such a protein^^ (Methanosarcina mazei CPD II 
photolyase) reported a time constant of —300 ms for this process 
at pH 8.0 {i.e., 200 OOOX slower than the 'slow' (-1.5 |is) direct 
protonation by D396(H) in AtCRYl). In insect CRYs, most of the 
FAD " appear to remain unprotonated for at least a couple of min- 
utes at pH 7.5^. These CRYs contain cysteine in the position corres- 
ponding to D396 in AtCRYl. As even the replacement of this 
cysteine by an aspartic acid^ did not lead to formation of FADH", it 
seems that the piCa of FADH" in insect CRYs is lower than 7.5. 

In a recent computational study in AtCRYl, Solov'yov and cow- 
orkers^^ predicted ultrafast protonation of FAD " by a two-step PT 
from TrpiH"+ via D396, the first step being PT from D396(H) to 
FAD ", and the second one PT from TrpiH"^ to D396". These PTs 
would become possible by structural rearrangements triggered by 
formation of the pair FAD " - TrpiH"^ that bring the groups 
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Figure 5 | Effects of ATP and pH on the absorption spectrum of FADqx in AfCRYl PHR. (a) Spectral change with increasing concentration of ATP 
(corrected for dilution due to the addition of ATP solution). Isosbestic points are situated at 317 nm, 362 nm and 41 1 nm. The inset shows ATP minus no 
ATP difference spectra scaled using 8445 = 11 300 M"^cm"^ for FAD ox- The sample contained 8 [iM protein, 50 mM phosphate buffer of pH 7.4 
and 0.5 MNaCl. (b) Comparison of spectra at pH 7.4 and pH 8.6 (50 |aM protein, 50 mM phosphate buffer of pH 7.4/50 mMTrisbuffer ofpH 8.6;both 
containing 0.5 M NaCl). Isosbestic points are situated at 315 nm, 364 nm and 409 nm. The inset shows the pH 7.4 minus pH 8.6 difference spectrum 
scaled in the same way as in panel (a). See Supplementary Fig. 3b for more details on the pH-dependence of the AfCRYl spectrum. 
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involved in the FT closer to each other. It was argued^^ that the 
observed 1.7 |is kinetics^^ was wrongly assigned to protonation of 
FAD " and may instead be attributed to deprotonation of D396(H) 
and protonation of Trp/ coupled to ET from Trp2H to Trpi" (the 
latter suggestion, however, is in contrast to the spectral features and 
amplitudes of the 1.5 |is kinetics observed in our study: such a pro- 
cess would yield virtually no absorption change at 403 nm where we 
observed a very prominent 1.5 jis phase). Nevertheless, analysis of 
our data (below) indicates that the ultrafast FT does indeed occur 
(contributing significantly in the absence of ATF), in competition 
with ET through the tryptophan cascade and subsequent microse- 
cond protonation of FAD " (strongly preferred in the presence of 
ATF). The products of the ultrafast FT process (FADH' - Trp/) and 
those of the classical ET through the Trp triad followed by micro- 
second FAD " protonation (FADH" - Trp3") recombine with time 
constants of 8.0 ms and 2.2 ms, respectively (see Fig. 4), unless Trp" 



radicals are reduced by tyrosine residues or by extrinsic reducing 
agents (see Supplementary Information, Section 4 for more details). 

The effect of pH on the stationary UV-Vis spectra of FADox and on 
the yield of photo-induced transients in the absence of ATF is a 
strong indication of a change in the dark state protonation of an 
amino acid (with pK^ ~ 7.4) in the immediate vicinity of the chro- 
mophore. We cannot exclude that the protonation states of more 
remote groups indirectly influence the properties of FAD, e.g., by 
affecting the solvent access to the FAD cavity. However, it appears 
much more plausible that the observed pH effects are essentially due 
to the protonation state of a residue directly interacting with the 
isoalloxazine ring of the FAD chromophore. Looking at the known 
crystal structure of AtCRYl FHR^°, possible candidates include 
aspartic acid residues D359, D390, D392 and D396. Because of its 
hydrophobic environment, the best candidate appears to be D396, 
the putative proton donor to FAD " The observed pK^ of —7.4 is 
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substantially higher than that of an aspartic acid in aqueous solution 
(piCa ^ 3.7). A hydrophobic environment destabilizes the negative 
charge of the deprotonated acid as compared to aqueous solution 
{i.e., increases the pK^ of the given acid). D359, D390 and D392 are 
exposed^^ to positively charged arginines (R360, R362 and R301, 
respectively) and are hence unlikely to have a pK^ higher than that 
of the free acid. D396 seems to be protonated (H-bond to the back- 
bone carbonyl oxygen of M381) in the crystal structure of AtCRYl 
PHR2° determined at pH 5.5, suggesting a pK^ > 5.5. An FTIR band 
attributed to light-induced deprotonation of D396 in AtCRYV^ and 
of its equivalent in the PHR domain of Chlamydomonas CRY^^ 
appears to be substantially more pronounced at pH 7.5^^ than at 
pH 8.0^^, which would be consistent with pK^ —7.4 obtained in 
our experiments. 

We propose that binding of ATP induces a further increase of pK^ 
of D396(H) from -7.4 to >9 (likely -10, as inferred from the 
decrease of signal amplitude at pH 9.6; Fig. 3), which assures that 
the protonated (uncharged) form of D396 prevails under physio- 
logical pH^^. An increase of the pK^ of D396 due to ATP binding 
could also explain the reported slower reoxidation of FADH" to 
FADox in the presence of ATP^^'^^, since FADH" needs to give off a 
proton in order to reoxidize back to FADox, which should be more 
difficult when D396 is protonated (see also Supplementary 
Information, Section 4). 

ATP and D396(H) lie — 10 A apart, separated by the isoalloxazine 
ring of FAD, thus the impact of ATP on the pK^ of D396(H) appears 
to be indirect. Possible mechanisms include (i) a decrease in polarity 
of the D396 environment due to ATP restriction of solvent access to 
FAD, (ii) an electrostatic effect from negative charges of the ATP- 
phosphates, and/or (iii) strains resulting from ATP stacking. The 
adenine moiety of ATP is sandwiched between L296 and Y402 near 
the adenine of FAD (Fig. 1). The 6-membered ring of the ATP 
adenine is pi-stacked with the phenol ring of Y402. D396, W400 
and Y402 are linked by helix otl6. In the crystal structures of 
AtCRYl PHR'° at pH 5.5, D396 appears to be protonated and H- 
bonded to the backbone carbonyl oxygen of M381 (al5), regardless 
of the presence of ATP. Since no crystal structures at higher pH are 
available, we can only speculate, consistent with the MD simula- 
tions^^, that Y402-ATP interaction likely puts constraints on the 
entire al6 helix, resulting in a modified environment for D396. 
For example, placing the al6 helix in a position where the D396 
carboxyl faces the electronegative backbone carbonyl oxygen and 
the sulphur atoms of M381 (al5) could increase the preference for 
protonation of D396 {i.e., its pK^) and the formation of an H-bond 
interconnecting the two helices al5 and a 16. 

We suggest two mutually non-exclusive mechanisms by which 
unprotonated (negatively charged) D396 might decrease the yield 
of long-lived radical pairs: 

(i) The ultrafast PT forming FADH' - Trpi" directly from FAD " - 
TrpiH"^ requires protonated D396. With unprotonated D396, 
this pathway is lost and ET through the tryptophan cascade is 
the only productive pathway competing with back ET in the 
pair FAD " - TrpiH +. 

(ii) By electrostatic interaction and/or by conformational changes, 
negatively charged (deprotonated) D396 might disfavour the 
first ET step (from TrpiH to * FADox) relative to the unpro- 
ductive decay of * FADox, and/or disfavour the second ET step 
(from Trp2H to TrpiH"^) relative to back ET in the pair FAD " 
- TrpiH"^, e.g., by repulsion between D396" of helix al6 (car- 
rying TrpiH) and the electronegative oxygen of the M381 back- 
bone carbonyl residing on helix al5 (carrying Trp2H), 
increasing the distance between TrpjH and Trp2H (weakening 
their electronic coupling). 

In addition to the prominent effects of pH and ATP on the ampli- 
tudes of the transient absorption signals discussed so far, there are 



also effects on the time course of the signals which are less pro- 
nounced, but significant. Thus, at 562 nm (see Fig. 2), addition of 
ATP at pH 7.4 increased the contributions of both the 200 ns phase 
(deprotonation of TrpsH'^ competing with unproductive back ET) 
and the 1.5 jis phase (establishment of the equilibrium FAD " + 
D396(H) ^ FADH" + D396") relative to the end amplitude (attrib- 
uted essentially to long-lived FADH" - Trp" pairs). In the framework 
of the suggested reaction scheme (Fig. 6), this implies that ATP 
favours the pathway involving ET through the tryptophan cascade 
over the ultrafast (<1 ns) PT pathway. Our observation that, in the 
absence of ATP, long-lived FADH" - Trp" can be formed with rather 
little TrpjH"^ deprotonation (200 ns) and 'slow' (1.5 |is) FAD " 
protonation supports the (co-) existence of the theoretically predicted 
ultrafast PT pathway^^. 

The effect of ATP on the branching between the two productive 
pathways may be qualitatively explained by the proposed increase of 
the piCa of D396 upon binding of ATP, which would decrease the 
driving force of PT from D396(H) to FAD " and hence disfavour the 
ultrafast PT pathway relative to ET through the tryptophan cascade. 
Of course, ATP may have other, more subtle effects on the branching 
ratio, such as conformational changes modifying ET and/or PT dis- 
tances and changes of driving forces and/or reorganization energies 
of ET. It is conceivable that the ATP-induced constraint to the al6- 
helix (see above) hinders the rearrangement of D396(H) required for 
the ultrafast PT. By hindering the solvent access to FAD, ATP might 
also affect other properties of FAD that are relevant to its ET and PT 
reactions. In this context, we considered the possibility that ATP 
binding might have a major effect on the reduction potential of 
FADox in AtCRYl but that does not seem to be the case (see 
Supplementary Information, Section 5). 

For a quantitative analysis and verification of consistency with the 
proposed reaction scheme (Fig. 6), we tried to globally fit all transient 
absorption data in Fig. 2 (i.e., signals at 403, 457 and 562 nm in the 
presence and absence of ATP at pH 7.4) by a fitting function based on 
the reaction scheme (see Methods and Supplementary Information, 
Section 7). It was assumed that an (ATP-dependent) fraction of the 
sample contained unprotonated D396 at the moment of excitation 
and could only undergo the reactions shown on white ground in 
Fig. 6, while another fraction contained protonated D396 and could 
undergo all the reactions shown. Reactions supposed to be faster than 
1 ns were not explicitly included in the model function but the pro- 
ducts of these reactions are contained in the early absorption changes 
of the model functions. Reverse ET through the tryptophan cascade 
was considered implicitly by introducing: (i) an effective back-ET 
rate constant for the pair FAD " - TrpsH"^ that includes thermal 
repopulation of the two preceding radical pairs, and (ii) an effective 
rate constant for formation of FADH" - Trpi" from FAD " - Trp3H"^ 
via repopulation of the first radical pair and PT assisted by D396(H). 
Absorption spectra of the relevant species (Fig. 2) were taken from 
the literature^^"^^ and from our own determinations (Supplementary 
Fig. 7). 

It turned out that the experimental data could be reasonably well 
fitted by our model (Supplementary Fig. 8b; fit parameters are sum- 
marized in Supplementary Tables 2 and 3). Of note, the intrinsic time 
constant of TrpjH"^ deprotonation (—350 ns) is slower than the 
apparent time constant from the multi- exponential fit (—200 ns) 
because the latter includes competing reactions (back ET and forma- 
tion of FADH" - Trpi" via repopulation of FAD " - TrpjH"^). 
Similarly, the intrinsic time constant of slow protonation of FAD " 
(—2.5 |is) is slower than the apparent time constant of — 1.5 |is that 
includes the reverse reaction. Further significant outcomes of the 
model-based fits are: 

(i) The fraction of the sample containing unprotonated D396 at 
pH 7.4 disappeared upon addition of ATP, in line with a strong 
increase of the pK^ of this aspartic acid. 
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(ii) In the fraction containing protonated D396, the branching ratio 
between the "ET through tryptophans" pathway (centre part of 
the proposed reaction scheme, Fig. 6) and the "ultrafast PT" 
pathway (right arrow in Fig. 6) changed from 2 : 1 in the absence 
of ATP to 7 : 1 in the presence of ATP. The "ultrafast PT" 
pathway is hence virtually negligible in the presence of ATP 
but it competes with the "ET through tryptophans" pathway 
in the absence of ATP. 

The fit results further indicate that the PT from D396(H) to FAD" " 
in the "ET through tryptophans" pathway is reversible with a for- 
ward-to-backward rate ratio in the order of 2, i.e., it has only a weak 
driving force. 

At first glance, it might seem puzzling that the same PT (from 
D396(H) to FAD ") can occur on such different time scales as <1 ns 
(before ET through the tryptophan cascade) and —2.5 |is (after this 
ET). There is, however, an important difference between the two 
situations: the nearby tryptophan W400 is in the radical cation state 
(TrpiH"^) in the former case, while it is neutral (TrpjH) in the latter 
case. By electrostatic interaction, the positive charge on W400 may 
strongly increase the otherwise weak driving force of PT from 
D396(H) to FAD ". Alternatively, the positive charge on W400 
may be specifically required for the ultrafast structural rearrange- 
ment that makes the PT possible. Once ET through the tryptophan 
cascade has occurred and W400 is neutral, the same structural 
rearrangement may be much slower. 

The modifications of the UV absorption spectrum of FADqx in 
AtCRYl upon decrease of the pH and/or upon addition of ATP 
resemble those observed upon binding of cyclobutane pyrimidine 
dimer (CPD) lesion to a DNA photolyase^^'^^. The latter effect has 
been explained by the strong electric dipole field of the CPD present 
in the substrate binding pocket^°. As adenine has a much weaker 
(~3.5X)^° electric dipole than the two nearly parallel thymines in 
the used CPD^^'^°, this explanation can hardly account for our obser- 
vations in A^CRYl. We rather suggest that the electric field of point 
charges (most plausibly the negative charge on D396 at pH > pK^ 
and possibly also the negative charges of ATP -phosphates) and the 
change of polarity of the flavin environment upon replacement of 
solvent in the binding pocket by adenine induce the spectral changes 
of FADox in AtCRYl. Such effects may also contribute to the spectral 
shifts accompanying substrate binding in DNA photolyase (see 
Supplementary Information, Section 2). 

Given the high affinity of AtCRYl for ATP (K^ < 2 |iM at pH 7.4, 
Supplementary Fig. 4) and the typical physiological concentration of 
ATP in plant cytoplasm (1.3 ± 0.3 mMy\ the ATP effects described 
here are unlikely to be used by the cells for regulation of crypto - 
chrome photoreceptor activity. ATP may rather be considered as a 
coenzyme of AtCRYl, the roles of which include increasing the yield 
of FADH" (regardless of possible pH fluctuations) and prolonging 
the lifetime^^'^^ of this putative signalling redox state of flavin (see also 
Supplementary Information, Section 4). A more plausible way to 
regulate cryptochrome light response in vivo is by modulation of 
the redox environment within the cell affecting the yield of FADH" 
(scavenging of Trp"/TyrO" radicals prevents their recombination 
with FADH"). 

We have shown that, in the presence of high concentrations of an 
extrinsic reductant, the in vitro quantum yield of long-lived FADH"/ 
FAD " can reach up to —14% (in the presence of ATP; the previous 
estimate^^ of only —2% had been obtained in the absence of ATP and 
at higher pH; see Supplementary Information, Section 8). This new 
value, combined with the prolonged lifetime of FADH" in the pres- 
ence of ATP^^'^^, seems much more congruent with the role of 
AtCRYl as photoreceptor and with formation of FADH" as the pro- 
cess responsible for cryptochrome signalling. 

Based on our analysis of the photoreduction pathways (Fig. 6) we 
would like to propose an alternative trigger of the signalling confor- 
mational change^^'^^ in AtCRYl. Unlike neutral FADH", the long- 



lived negative charge of D396" should induce a strong perturbation 
of the hydrophobic protein environment and is hence much more 
likely to trigger a major conformational change (leading to transduc- 
tion of the initial light signal), e.g., by electrostatic repulsion between 
D396" (helix al6) and the electronegative backbone carbonyl oxy- 
gen of M381 (helix al5), to which protonated D396 is hydrogen- 
bonded prior to photoexcitation^°. Long-lived D396" is formed only 
by the ATP-enhanced "ET through tryptophans" pathway (centre of 
the reaction scheme. Fig. 6) but not by the "PT first" pathway, in 
which D396" is promptly protonated by TrpjH"^. Another effect of 
ATP, namely the up-shift of the pK^ of D396 from 7.4 to >9, is 
crucial to make this photoswitch efficient as it assures that virtually 
all cryptochromes are turned "off (D396(H) state) in the dark under 
physiological conditions. Future experiments and computer simula- 
tions should test these hypotheses. 

Methods 

Reagents and materials. Water was de-ionized (> 18.2 MQcm specific resistance) by 
Elga ultrapure water purification system. Adenosine 5 '-triphosphate disodium salt 
(>99%), Imidazole (>99%) and [Ru(bpy)3]Cl2 (99.95%) were purchased from Sigma 
Aldrich Chemicals. Tris buffer was prepared by titration of 
tris(hydroxymethyl)aminomethane (EUROMEDEX, molecular biology grade) by 
HCl (Carlo Erba Reagenti,p.a. grade). Phosphate buffer was prepared from Na2HP04 
ip.a.) from Merck and NaH2P04 2H2O (>99%) from Prolabo. Triton X-100 was 
purchased from Prolabo. Quartz cells for spectroscopic experiments with self- 
masking solid black walls and sample chambers of 2 X 8 X 10 mm (width X height X 
length) with 4 clear windows were obtained from Starna. 

Protein preparation. The plasmid coding full-length AfCRYl with a His Tag was 
inserted into the BacPAK™ 6 baculovirus expression system (Clontech Laboratories, 
Inc.). Sf21 insect cells were inoculated with the baculovirus and incubated at 28°C. 
After 48 h, the cells were spun down (2 600 X g) and lysed in a lysis buffer (50 mM 
Tris buffer at pH 7.6, 0.5 M NaCl, 1% Triton X-100) for 60 min on ice (lysis buffer 
volume ~ lOX the volume of the cell sediment). The cell lysate was spun down ( 19 500 
X g) and the supernatant was incubated with NTA resin from Qiagen (0.5 to 1 mL per 
50 mL) at 4°C for 30 minutes. The resin was washed with 50 mL of the wash buffer A 
(50 mM Tris buffer at pH 7.6, 0.5 M NaCl, 10 mM imidazole), with 10 mL of the 
wash buffer B (50 mM Tris buffer at pH 7.6, 0.5 M NaCl, 20 mM imidazole) and 
eluted with elution buffer (50 mM Tris buffer at pH 7.6, 0.5 M NaCl, 350 mM 
imidazole). Imidazole was removed on Micro Bio-Spin (Bio-Gel P-6) size-exclusion 
chromatography columns from Bio-Rad. 

Gene coding the AtCRYl lacking the CCT domain was inserted into pGEX-6P-l. 
The protein was expressed in E. coli and purified with glutathione Sepharose followed 
by Mono Q and Sephacryl S-300 after tag cleavage. The protein was dissolved in a 
buffer solution containing 50 mM Tris-HCl (pH 8.0), 50 mM NaCl, and 5% (vol/vol) 
glycerol and deep-frozen ( — 80°C). Before use, the buffer was exchanged to 50 mM 
phosphate (or Tris) buffer of desired pH (and 0.5 M NaCl) using Micro Bio-Spin 
(Bio-Gel P-6) size-exclusion chromatography columns. 

UV-Vis and transient absorption spectroscopy. UV- Vis spectra of protein solutions 
and the reference compound [Ru(bpy)3]Cl2 were recorded using the Uvikon XS 
spectrophotometer (Secomam). Protein spectra were recorded prior to, several times 
during and after each experiment in order to verify there was no significant 
accumulation of FADH' or deterioration of the samples by multiple laser flashes. 

In all transient absorption experiments, the samples were excited at 355 nm by a 
Nd : YAG laser (Continuum Leopard SS-10, pulse duration of 100 ps, repetition rate 
0.25 Hz). The vertically polarized laser beam was passed through attenuation filters 
and shaped by lenses so that the cross section of the beam covered the whole 8 X 
10 mm cell window. Laser energies (as a control for the actinometric quantum yield 
calculation) were measured behind a cell filled with H2O using an energy meter 
(Gentec QE25SP-H-MB-D0). 

Experiments on timescales <20 |j,s were performed using a modification of the 
setup described previously^'*. Monitoring light was provided by the following light 
sources: 

403 nm - 105 mW Laser diode SHARP GH04P21A2GE, 

457 nm - 50 mW DPSS cw- Laser Cobolt Twist, 

562 nm - 25 mW DPSS cw-Laser Oxxius 561-25-COL-002. 

The monitoring light beam was perpendicular to the excitation beam and shaped 
such that the central part of the beam passed through the 2X8 mm cell window 
along the 10 mm path. In order to minimize sample excitation by monitoring light, 
attenuation filters and a rotating blade with a small opening was introduced between 
the light source and the sample, chopping the monitoring light to rectangular pulses 
of 140 |is duration (separated by —30 ms dark intervals). The energy of one such 
140 |is pulse of measuring light received by the sample was in the order of 1 [l] (beam 
cross-sectional area of ~5 mm^), which is negligible compared to ~4 mj/cm^ of the 
excitation flash. 

Nevertheless, the repeated monitoring light pulses could lead to some loss of dark 
state (FADox), particularly for monitoring light at 457 nm, where FADqx absorbs 
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most strongly (s = 9.5 mM~^cm~^). Taking into account an AtCRYl concentration 
of —45 |iM, a maximal quantum yield of — 15% (in the presence of ATP) of radical 
pairs at 10 [is after excitation (Supplementary Information, Table 4) and recom- 
bination of —85% of these radical pairs during the dark time of 30 ms between two 
monitoring light pulses (Fig. 4), we estimate that at most 6.5% of the dark state could 
be lost due to the monitoring light during accumulation of one package of 16 transient 
absorption signals (total acquisition time of — 1 min). In fact, the observed accu- 
mulated loss of dark state (as estimated from the transmission increase at 457 nm) 
due to the combined effects of monitoring light and excitation flashes never exceeded 
2% indicating that further radical pair recombination occurs in the time window from 
30 ms to 1 min. 

The system was synchronized so that the excitation pulse came 30 |is after the rise 
of the monitoring light. A X/2 plate (Melles Griot 99E02 WRA 0000) and a polarizer 
(COLPA 10.0-425-675 set to 54.7° with respect to the vertical direction) were inserted 
between the light source and the sample to turn the polarization of the monitoring 
light source so that the excited molecules were probed at magic angle to the excitation 
light polarization (elimination of anisotropic artefacts). The position of the cell with 
respect to the excitation beam was unchanged for all experiments in order to guar- 
antee comparable excitation conditions. 

The detection system for the monitoring light consisted of a Si photodiode 
(Alphalas UPD-500-UP, rise time 500 ps, spectral range 170-1 100 nm, sensitive area 
1 mm^), connected via an electronic signal amplifier (Femto HCA; DC-325 MHz, 
28 dB) to a digital oscilloscope (Tektronix DSA602A, with 11A52 plugin set to DC- 
20 MHz). Interference filters with transmittance maxima at the wavelengths of the 
monitoring light were placed in front of the photodiode to block stray light and 
possible sample fluorescence induced by the excitation pulse. The oscilloscope was 
triggered by a small fraction of the excitation pulse through a fast photodiode. The 
kinetic curves are averages of 64 to 128 signals recorded in packages of 16 signals with 
several minutes dark interval between the packages. 

For experiments on time scales >20 |j,s, a 100 W-tungsten-halogen lamp was used 
as a light source, combined with interference filters with transmission maxima at 410, 
450 or 560 nm and spectral bandwidths of —5 nm placed between the lamp and the 
sample. The probe light intensity and its changes associated with the kinetics of the 
light-induced processes were monitored behind the sample by a Si photodiode (FND 
lOOQ from EG&G, protected by another interference filter transmitting at the 
respective wavelength) connected to an amplifier (Tektronics AM502, upper band- 
width limit set to 1 or 10 kHz) and the amplified signal was fed to a digital oscilloscope 
(Tektronix DSA602A, with 1 1A33 plugin). The kinetic curves are averages of 4 signals 
(with 4 s intervals between two excitation pulses allowing for sample recovery 
between the individual pulses). 

The samples were air- saturated and constantly kept on ice or in cooled sample 
holders (keeping the temperature between 0 and 2°C) to minimize the risk of protein 
aggregation and degradation. The pH values of the buffers refer to this temperature 
range. 

Quantum yield determination. Quantum yields of the transient radical pairs were 
determined using the [Ru(bpy)3]^^ actinometer^^-^^-^*^. 44.7 [iM [Ru(bpy)3]Cl2 (in a 
cell identical to that containing AtCRYl) was excited under the same conditions as 
the protein in the experiment, the results of which are shown in Fig. 2. The average 
signal amplitude of this actinometer allowed estimating the incident photon fluence 
per excitation pulse and the concentration of excited FADqx in the AtCRYl samples. 
Quantum yields of the different radical pairs in AtCRYl were calculated by referring 
their respective concentrations (obtained from the analysis of the transient 
absorption data) to the concentration of excited FADqx (see Supplementary 
Information, Section 8). 

Spectroelectrochemistry. Spectroelectrochemical titrations of AtCRYl PHR were 
performed essentially as described in Ref 47. See Supplementary Information, 
Section 5. 

Fluorescence spectroscopy. Fluorescence measurements (Supplementary 
Information, Section 1) were conducted using the FluoroMax®-4 spectrofluorometer 
from Horiba Scientific. Excitation was at 445 nm (bandpass 10 nm), detection was at 
540 nm (bandpass 5 nm). 

Global analysis of the transient absorption kinetics. For a phenomenological 
description, the six experimental curves in Fig. 2 were globally fitted by multi- 
exponential decays. A satisfactory fit (see Supplementary Fig. 8a) was obtained with 
the tri-exponential function: 



AAA,+ (t)=AAu,±e-^^' + AA2, 



-k2t 



+ AA3. 



-kit 



(1) 



where AA^ + denotes absorption change at wavelength /I in the presence ( + ) or 
absence (-) of ATP, t is time, AA^ to AA3^^ + are the amplitudes of the individual 
kinetic phases at the given wavelengths and ATP condition, and ki to the 
exponential decay constants common to all six signals. The constants AAo^x,± account 
for absorption changes that are much longer-lived than the 10 |as time window of the 
fit. The best fit yielded ki = 4.1 X 10' s'K k2 = 5.3 X 10' s"^ and k^ = 6.6 X 10=^ s"^ 
with amplitudes listed in Supplementary Table 1. Allowing the decay constants k^ to 
depend on the ATP condition did not significantly improve the quality of the fit. The 
k2 and the k^ components were in first approximation assigned to deprotonation of a 
tryptophan cation radical TrpH'"^ and protonation of FAD ", respectively (see 



Results). The ki component is close to the time resolution limit of the set-up and 
reflects the instrument response to processes that are too fast to be correctly 
monitored in the present experiments. The amplitudes AA^ ^ + were not used for 
further assignments. 

The six experimental curves in Fig. 2 were also fitted by a function (equation (2)) 
that is based on the reaction model shown in Fig. 6. For the sake of simplicity, we 
introduce abbreviations for the states that could be observed at the present time 
resolution: 

A = FAD " - Trp3H +; B = FAD " - Trp3'; C = FADH" - Trp3'; C = FADH" - 
Trpi' and D = FADqx - TrpH (ground state). Rate constants /cab> ^ac'> ^ad> ^bc and 
kcB refer to reactions A^B, A^C',A^D, B^C and C ^ B, respectively. The 
species kinetics observable on the time scale of our experiments (—10 ns-10 |is) can 
be described as a superposition of three contributions corresponding to three dif- 
ferent partial pathways (named a, P and y in Supplementary Chart 1): 

• Proteins that contained unprotonated D396 at the moment of excitation and 
escape direct back FT in the first radical pair FAD " - TrpiH'"^ by FT through 
the tryptophan cascade. Protonation of FAD " cannot occur. The initial (on our 
time scale) concentration of the pair FAD " - Trp3H (A) in the proteins fol- 
lowing this pathway is denoted c^. 

• Proteins that contained protonated D396 at the moment of excitation and 
undergo FT through the tryptophan cascade (rather than immediate back FT 
or D396(H) -assisted ultra-fast PT in the first radical pair FAD " - TrpiH +). The 
initial (on our time scale) concentration of the pair FAD " - Trp3H'"^ (A) in the 
proteins following this pathway is denoted Cp. 

• Proteins that contained protonated D396 at the moment of excitation and 
undergo D396(H) -assisted ultrafast PT in the first radical pair FAD " - 
TrpiH'"^ (rather than direct back FT or FT through the tryptophan cascade). 
The initial (on our time scale) concentration of the pair FADH' - Trp/ (C) in the 
proteins following this pathway is denoted Cy. 

The differential equations governing the species kinetics, their analytical solutions 
and the conversion to absorption change kinetics are given in the Supplementary 
Chart 1. A term AA; e~'^^ ^ was added to account for absorption changes that are too 
fast to be correctly resolved in the present experiments (see explanation to equation 
(1) above). The following equation was obtained: 

AAx,±{t)=AA',^^e-''^' 

c.,±|A£A,ie-(<±-'^-'^±)^ + A£B,;i 

-kAB,± ^-(k[^-k,,,^)t ^AB,± 



+ 



+ cp, + |A£A,;e"''2^±' + A8B,x 

^AB,± (/C2 + -/C3 + +/CbC,: 



^ kAB,±kBC,± ^-k[^t I 



-kAC,± 



(2) 



^3,+ (^2,+ ^3,+ ) 



+ Ae, 



^AB,+ ^BC,± 



^AC',±\ - 



.^2,+ (^2,+ ~^3,±) ^2,± 

^AB,+ ^BC,± Q-l^i+t ^ ^AB,+ ^BC,± _^ ^AC',± \ I 
3,± (^2,+ ~^3,±) ^2,+ ^3,± ^2,+ J J 



+ Cy, + Aec,. 



The indices + and — refer to the presence and absence of ATP, respectively. k2 



^AB,± + ^AC,± + ^AD,± and/c3 + = /cbc,± + kc 



Aec 



- ^¥ADo^,l + ^TrpH + -a ' 



-eFADox,A + STrp ,; and 



ADH A ~^FADox,A + eTrp A' ^TrpH was Omitted because neutral tryptophan in 
the ground state does not absorb at the wavelengths studied. C has the same 
absorption as C. 

In globally fitting the six experimental curves in Fig. 2 with equation (2), ki was 
fixed at the value obtained in the fit with equation (1). All other rate constants were 
free and allowed to depend on the ATP condition. Absorption coefficients were either 
fixed at the values from the spectra shown in Fig. 2 (for FADqx, TrpH'"^ and Trp' at all 
wavelengths and for FADH' at 457 and 562 nm) or allowed to vary in a limited range. 
FAD " spectra in Ref 35, as well as in other studies on insect cryptochromes (where 
FAD " is much longer lived than in plant cryptochromes), have apparently been 
obtained by adding up- scaled photoreduction difference spectra to the spectrum of 
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FADox- Because of an uncertainty with respect to the exact shape and scaling of these 
spectra, we have decided to use directly the photoreduction difference spectra of 
Drosophila cryptochrome from Ref. 39 to determine the ratio of [^p^o a ~ ^faDox,^] 
the three wavelengths (403 nm/457 nm/562 nm): 15/-12/1 (note that because the 
blue absorption band of FADqx in AtCRYl {Imax — 445 nm) is blue-shifted by 5 nm 
with respect to insect CRY {Imax = 450 nm), AA462 instead of AA457 was read out 
from the figure). Maintaining this ratio was a condition in our fit and the absolute 
scaling at 403 nm was a free parameter within the range given by published spectra of 
protein-bound FADox and FAD " 3539.48,49. 5qqq M'^cm-^ < 
[%AD - 403 nm ~ ^FADox,403 nm] — 7000 M"^ cm"\ The dip close to 400 nm in our 
FADH" spectrum (Fig. 2 and Supplementary Fig. 7) is less pronounced than in FADH' 
spectra from other flavoproteins^°. We cannot exclude that this is due to some con- 
tamination of our spectrum by FADH" or by a species that had been oxidized by Trp" 
during photoreduction of FAD ox- Therefore, Sp^DH 403 nm ~ ^faDox,403 nm was allowed 
to vary between -1500 M"^cm"^ (value in Fig. 2) and -3000 M"^cm"\ All fits were 
done with user- defined functions in the Nonlinear Curve Fit algorithm of Origin 
(version 8.6). 
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